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Strong repulsive interactions in a one-dimensional electron system suppress the exchange coupling 
J of electron spins to a value much smaller than the Fermi energy Ep. The conventional theoretical 
description of such systems based on the bosonization approach and the concept of Tomonaga- 
Luttinger liquid is applicable only at energies below J. In this paper we develop a theoretical 
approach valid at all energies below the Fermi energy, including a broad range of energies between 
J and Ep. The method involves bosonization of the charge degrees of freedom, while the spin 
excitations are treated exactly. We use this technique to calculate the spectral functions of strongly 
interacting electron systems at energies in the range J <^ e <C Ep. We show that in addition to the 
expected features at the wavevector k near the Fermi point kp, the spectral function has a strong 
peak centered at = 0. Our theory also provides analytical description of the spectral function 
singularities near 3kp (the "shadow band" features). 

PACS numbers: 71.10.Pm 



I. INTRODUCTION 

Recent experiments with quantum 

wires^>^>^>^>^>^>^>^i^>^ and carbon nanotubes^^'^^ have 
stimulated theoretical interest in transport properties of 
one-dimensional systems of interacting electrons. It is 
now widely accepted that in one dimension interacting 
electrons form the so-called Luttinger liquidJ^ii^ The 
main signature of the Luttinger liquid — the power-law 
behavior of the tunneling density of states — has recently 
been observed in experimentsj ^^'^^i^^ Another well- 
known predictioni^iiiiiS. of the Luttinger-liquid theory 
is that the conductance of a quantum wire connecting 
two-dimensional leads should be quantized in units of 
2e^//i, regardless of the interaction strength. Although 
the quantization of conductance is routinely observed 
in modern experiments, careful recent measurements 
show significant deviations^iii^iii^i^ii^ from perfect 
quantization in the regime of very low electron density, 
where the effective electron-electron interactions are 
very strong. 

The applicability of the Luttinger-liquid theory is not 
expected to be limited to weak interactions. On the other 
hand, the properties of the system do change significantly 
when the interactions become strong. It is well known 
that at low energies one-dimensional electron systems 
support separate charge and spin excitation modes prop- 
agating at different velocities'^ Vp and Va- Accordingly, 
the Luttinger-liquid theory describes the low-energy ex- 
citations of the system by two bosonic fields with lin- 
ear dispersion, propagating at velocities Vp and Va-- The 
applicability of such a description is limited to energies 
small compared to the bandwidths of the charge and spin 
excitations Dp a- TiVp_„n, where n is the electron den- 
sity. In the non-interacting case, both velocities coincide 
with the Fermi velocity vp, so Dp — D^r ^ Ep. In the 
presence of weak repulsive interactions the velocities are 



renormalized, so that v^- < Vp, but both velocities re- 
main of order vp ■ At strong interactions the spin mode 
velocity Va is strongly suppressed, Va- ^ Vp. In this case 
D„ ^ Ep < Dp, and the Luttinger-Hquid theory is appli- 
cable only to phenomena in which all the relevant energy 
scales are smaller than D^- 

A number of recent theory 

papers22^^i23^>25^i2i^^ addressed the physics 
of strongly interacting electrons beyond the range of 
applicability of the Luttinger-liquid theory. Penc et 
^^ ^20,21,22 gt^fjied the tunneling density of states and 
spectral functions of the one-dimensional Hubbard 
model at energies in the range D^ <C £ <C Dp and zero 
temperature. Cheianov and Zvonarev^^'^"* and Fiete and 
Balents^^ explored the so-called spin incoherent regime^ 
D„ <C T ^ Dp and found an enhancement of the tun- 
neling density of states at energies e > T. Conductance 
of the quantum wire entering the spin-incoherent regime 
was predicted^liS^ to show behavior similar to the 
anomalies observed in experimentSi^i^i^iii^i^iiS, 

Despite the recent theoretical successes in treating 
strongly interacting one-dimensional electrons, at present 
there is no regular theoretical technique that can be ap- 
plied to a broad class of problems and is not limited 
to the exactly solvable models i^i^ii^iS^iiH In particular, 
the bosonization technique commonly used to justify the 
Luttinger-liquid pictur o^'^'^^ is applicable only at energies 
below the spin bandwidth Da-. In this paper we general- 
ize the bosonization technique to all energies below the 
charge bandwidth Dp. Our method treats the spin ex- 
citations exactly, but applies bosonization to the charge 
excitations. It is thus applicable to a broad class of strong 
interactions, and is not limited to short-range coupling 
required for the existence of the exact solutions. 

We apply our technique to the calculation of the spec- 
tral functions and the tunneling density of states of 
strongly interacting one-dimensional electron systems. 
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Unlike many of the earlier treatmentS )^°i^^'^^i^^'^ our 
calculations can be applied to systems with long-range 
interactions, such as quantum wires. In addition, our 
theory is valid in a broad range of temperatures: we ob- 
tain the spectral functions at zero temperature, in the 
spin-incoherent case T ^ D^, and also interpolate be- 
tween these regimes. Furthermore, our approach pro- 
vides a clear physical picture of the enhancemeni>2£i2^i^ 
of the tunneling density of states v{£) at |e| ^ Ep- 

Our approach is introduced in Sec. [TTl where we derive 
the expression for the electron field operators at energies 
below Dp by bosonizing the charge modes while treating 
the spin excitations accurately, as excitations of an ef- 
fective Heisenberg spin chain with a small exchange con- 
stant J. At energies below Dc J the spin excitations 
can also be bosonized. This is accomplished in Sec. IIIIl 
where we also demonstrate that the standard bosoniza- 
tion expression for the electron operator a^^'^^ is recovered 
in our approach when all important energy scales are be- 
low J. Calculation of various physical properties of the 
system requires knowledge of the electronic Green's func- 
tions discussed in Sec. lIVI In the most interesting regime 
of energies |e| ^ J the Green's functions are expressed in 
terms of certain equal-time correlation functions {q) of 
the Heisenberg spin chain. Their behavior is important 
for understanding the electronic transport at |e| ^ J; it 
is discussed in Sec. IIVI Wc calculate the spectral func- 
tions at energies J ^ |e| ^ Ep in Sec. fVl and show that 
their dependence on the wavevector k has a Gaussian 
peak centered at fc = 0, which determines the behavior 
of the tunneling density of states. In addition, we find 
power-law singularities in the spectral functions at k near 
fci?, as well as the shadow band features^^ near ikp. We 
conclude the paper with the discussion of our results in 
Sec. IVIl A brief summary of some of our results was 
reported in Ref. [3l|. 

II. BOSONIZATION OF CHARGE 
EXCITATIONS 

The most experimentally relevant one-dimensional sys- 
tem of strongly interacting electrons is realized in GaAs 
quantum wires i^i^i^iii^iiSiiiS In these device the spectrum 
of electrons is quadratic, and the system can be described 
by the standard Hamiltonian 

f 

H = j i)\{x)dlij-f{x) dx 

+ IJJ ijl{x)^l{y)V{x - y)il;p{y)^^{x)dxdy. (1) 

Here rrig is the effective mass of electrons, ip'y{x) is the 
annihilation operator of electron with spin 7 (summation 
over repeating spin indices is implied), potential V(x~y) 
describes the interaction between electrons, and dx = 
dj dx. 

The assumption of quadratic spectrum in the Hamil- 
tonian (1) is introduced for simplicity, and most of the 



results we obtain apply to a generic spectrum. In particu- 
lar, our theory is applicable to the one-dimensional Hub- 
bard model, with an important exception of the half-filled 
case, where the charge excitation spectrum is gapped. 

The effect of interactions on the low-energy proper- 
ties of the system is quantified by the parameters r]f — 
V{^)l'hvp and r\h = V{2kp)/hvp, where V{q) is the 
Fourier transform of the interaction potential. Parame- 
ter rjf controls the amplitude of forward scattering of two 
electrons at the Fermi surface. Positive value of ry/ leads 
to the enhancement of the velocity of charge excitations 
Vp over the Fermi velocity. Parameter rj}, controls the am- 
plitude of backward scattering of two electrons. Strong 
backscattering impedes propagation of spin excitations 
through the system and leads to the suppression of the 
spin velocity v„. In the case of short-range interactions, 
the two parameters are of the same order of magnitude. 
On the other hand, in quantum wires the electrons inter- 
act via long-range Coulomb repulsion, which is usually 
screened at a large distance d by a metal gate. In this 
case rjf/rih ~ In(nd). Throughout this paper we assume 
strong backscattering, rjb^ 1- 

A. Effective Hamiltonian 

In the limit r^f, — s- cx) collisions of two electrons with 
opposite spins result in complete backscattering. In this 
case the processes of spin exchange are completely sup- 
pressed, and the energy of the system no longer depends 
on the spin degrees of freedom. To find the energy of any 
state, one can assume that all the spins 7 —'\, or, equiv- 
alently, assume that the fermions in the Hamiltonian ([T]) 
are spinless, ■07 (^) ~^ ^(a;). The resulting Hamiltonian 

Hp - -^J^\x)dl^{x)dx 

+ 'i\x)-^\y)V{x - y)-^[y)-^{x) dxdy (2) 

describes the charge excitations in the system. Each 
eigenstate of the Hamiltonian ^ for a system of N elec- 
trons is a degenerate multiplet of 2^ spin states. 

In this paper the spinless fermions ^(x) will be referred 
to as holons. By construction their number equals the 
total number of electrons, 

^\x)^{x) ^^{x)iI^-^{x)+iIj\{x)iIj^{x). (3) 

In the limit of short-range coupling, V{x — y) = Vq5{x — 
y) , the interactions in the Hamiltonian ([2]) disappear due 
to the Pauli principle, ^{x)'^{x) = 0, and holons become 
free fermions. This fact is well known in the theory of 
the Hubbard model at strong interactionsi^ 

The physical picture of one-dimensional electrons in 
the limit of strong backscattering is illustrated in Fig. [T] 
Due to the strong repulsive interactions, electrons cannot 
pass through each other. As a result electron I is always 
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FIG. 1: One-dimensional electrons at strong backscattering. 
Each electron is confined to the region of space between the 
neighboring particles. Shaded regions represent strong repul- 
sive interactions. 

confined between electrons / — 1 and l + l, and their spins 
cannot move through the system. 

At strong but finite repulsion the amphtude of forward 
scattering of electrons does not vanish. Such processes 
give rise to a weak exchange of the spins at the neighbor- 
ing sites, e.g., I and Z -I- 1. To leading order the couphng 
of the next-nearest neighbors can be neglected. The sym- 
metry with respect to spin rotations dictates the form of 
coupling between the spins: 

H, = Y,JSi- Si+i. (4) 
I 

The exchange constant J is positive, as otherwise the 
ground state of a system of one-dimensional electrons 
would have been spin-polarized, in violation of the Lieb- 
Mattis theoremi^ 

Thus we conclude that the Hamiltonian of a one- 
dimensional system of strongly interacting electrons can 
be written as a sum Hp + Ha. In the case of the 
Hubbard model this was first noticed by Ogata and 
Shiba,'^^ who discovered that in the limit of strong re- 
pulsion the Bethe ansatz ground state of the system fac- 
torizes into a direct product of the ground state of non- 
interacting fermions (holons) and the ground state of the 
Heisenberg spin chain ([¥]). In the context of quantum 
wires the description based on the Hamiltonian Hp + H^ 
was used in Refs. [2711281 . The long-range nature of the 
Coulomb repulsion between electrons in a quantum wire 
results in the exponential suppression of the exchange 
constant^^>^ 

J = J*expf-4^). (5) 

Here the prefactor J* ^ EpinaB)"^^^ , the Bohr radius 
is defined as as — e?i^/mee^, and e is the dielectric con- 
stant. 

It is worth mentioning that one can replace the Hamil- 
tonian U]) with the sum of two independent Hamiltonians 
and ([3]) only at sufficiently low energies. Indeed, the 
backscattering amplitude for two electrons with wavevec- 
tors ±fc is V{2k)/hvk, and tends to zero at k oo for 
any reasonable interaction potential. (Here Vk — hk/m 
is the velocity of an electron with wavevector fc.) Thus at 
high energies the initial assumption of strong backscat- 
tering is violated. To find the region of applicability of 
our low-energy theory, one can estimate the correction 
to the exchange constant (O caused by the fact that 
at higher energies the distances between electrons fluc- 
tuate, and the density n is no longer constant. Given 



that the rigidity of the Wigner crystal is due to the 
Coulomb repulsion between the electrons, we estimate 
Sn ~ y^m^naBlel/h, and the exchange J acquires sig- 
nificant corrections in the presence of excitations with 
energies \e\ > [Tm)^ /rrie ~ Ep- Thus our subsequent re- 
sults are valid up to energies of order Ep, rather than the 
somewhat higher energy scale Dp ^ Ep y^hi{nd)/naB- In 
the case of short-range interactions the scales Dp and Ep 
are of the same order of magnitude. 

B. Electron creation and destruction operators 

The effective Hamiltonian ([2]), (|4]) is defined in terms 
of the holon field operators '^(x) and the spin operators 
Si, rather than the original electron operators V'T,i(^)- 
In order to apply the effective theory (21), (HI) to prob- 
lems formulated in terms of electrons (e.g., calculation of 
the spectral functions) we need to establish the relations 
between the electron operators and the new variables. 

In the context of the Hubbard model this issue was 
addressed by Penc et al.^^ who used the definition 

V't(0) = *t(0)4^. (6) 

Given the relation ^ between the densities of electrons 
and holons, creation of an electron at point must be 
accompanied by creation of a holon. In addition, when a 
new particle is added to the system of N electrons, the 
spin chain ^ acquires an additional site. This is ac- 
counted for by the operator Zq ^. By definition, operator 

zj ^ adds a new site with spin 7 to the spin chain between 
the sites I — 1 and I. 

Despite the fact that the rule (O leads to a number of 
correct results when applied carefuUy^^i^ii^ it is not a 
completely satisfactory expression of an electron creation 
operator in terms of the charge and spin degrees of free- 
dom. In particular, the generalization of Eq. ^ to x ^ 
is not straightforward.^^ The origin of the difficulty lies 
in the fact that unlike the electrons and holons, the spins 
in the Hamiltonian ^ are not assigned to specific points 
in space. When an electron is created at point x, the 
additional site in the spin chain ^ appears at I = l{x), 
where 

l{x)= r^\\v)^{y)dy (7) 
J —00 

is the number of electrons (or holons) between —00 and 
point X. We shall therefore define the electron creation 
and annihilation operators as 

^\{x) = Z{^^^^^\x), (8a) 
ip^{x) = ■^{x)Zn^^)^^. (8b) 

The most important difference between the expressions 
^ and ([5a|) is that the latter explicitly accounts for the 
fact that the spins are attached to electrons. Thus despite 
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the apparent separation of the charge and spin degrees 
of freedom in the effective Hamiltonian (H]), (51), the elec- 
tron creation operator (|8ap does not factorize into a prod- 
uct of two operators acting on only charge or only spin 
variables. "^^ In Appendix \X\ we show that our operators 
(IH) satisfy the appropriate anticommutation relations. 



C. Bosonization of holon operators 

In this paper we are interested in the properties of 
strongly-interacting one-dimensional systems at energies 
well below the Fermi energy, \e\ <^ Ep. In this case the 
dynamics of the charge degrees of freedom described by 
the Hamiltonian ([2]) simplifies dramatically. Indeed, it is 
well known that the low-energy properties of a system of 
interacting spinless fermions are accurately described by 
the Tomonaga-Luttinger modelfi^ii^ 



+ K-\d.,(^f] dx. 



(9) 



Here (}> and are bosonic fields satisfying the com- 
mutation relations [(j){x),dy9[y)\ — TTi6{x — y). Field 
6 is related to the momentum density of the system, 
p{x) — hndxO{x), whereas </> is defined in terms of the 
density of fermions, 
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(10) 



Here we have introduced the Fermi momentum of the 
holons kp. Considering that the total electron density 
equals the holon density, Eq. ([3]), we have n = 2kp /i: = 
k^/iT. Thus the Fermi wavevectors of holons and elec- 
trons are related by 



2k 



F- 



(11) 



IS 



Finally, the parameter K < 1 in the Hamiltonian 
the so-called Luttinger- liquid constant. 

Transformation of the Hamiltonian Hp from the 
fermionic form ([2]) to the bosonic form ([9]) is accom- 
plished via bosonization procedure /^'^^ in which the 
fermion operators ^E* are expressed in terms of the bosonic 
fields (f) and 6. At the first step one notices that at low 
energies the properties of one-dimensional Fermi systems 
are dominated by excitations near the two Fermi points. 
Particles near each of the Fermi points propagate in one 
of two possible directions, right or left. Thus the fermion 
operator is presented as a sum of two chiral fermions, 



(12) 



where the operators r and l destroy fermions with 
wavevectors near the right and left Fermi points, respec- 
tively. 

The chiral fermion operators are bosonized following 
the prescription 



*i?,L(a;) 



1 



(13) 



where a is a short-distance cut-off. 

Using the relations ([7]), pU)). (fT^ . and (fT^ . we express 
the electron annihilation operator (|8bp in the form 



■07 (2;) 



-ie(x) 



V27ra 



^^i[k%x+,p{x)] _|_ g-i[fc^x+0(x)] 



i=i;[k'^x+^(x)] 



(14) 



Unlike the original formula (jSbp . this expression is valid 
only at low energies, |e| <C Ep- Its advantage is that 
the charge modes are now presented in the form of non- 
interacting bosons (O and can be treated rather easily. 
We apply the expression (fT4|) to the calculation of the 
electron Green's functions at J, T ^ lei ^ Ep in Sec. lIVI 



III. BOSONIZATION OF SPIN EXCITATIONS 

As we discussed in Sec. HI the conventional Tomonaga- 
Luttinger theoryi^ii^ of low-energy properties of one- 
dimensional electron systems is based on the idea of 
bosonization of electron operators. Mathematically this 
corresponds to applying the procedure (|13p to annihila- 
tion operators of electrons with spins 7 



1pLi{x) 



e'^""' ^[^^(x)-e,{x)] 



^/2^ 



^^±^[4,a(x)'-e,(x)] 



\/27rQ! 



(15a) 



(15b) 



Here (pp^^ and 9pu are the bosonic fields describing the 
charge and spin excitation modes of the system. In terms 
of these fields the Hamiltonian of interacting electrons ^ 
takes its bosonized low-energy form Hp + H^ with 



H„ = 



flVn 



27r 



^ / [Kpid^Bpf + K;\d.x(l)pf] dx, (16a) 



^ j [KAdx0.f+K-\d,<j,,f]dx 



251 



(27ra)2 



cos 



\/8 0<,(a;) 



dx. 



(16b) 



Here the matrix element gi±^ accounts for the processes of 
backscattering of two electrons with opposite spins. The 
respective sine- Gordon term in the Hamiltonian (|16b[) is 
marginally irrelevant. In the absence of magnetic field, 
the SU(2) symmetry of the problem requires^^ that when 
gi^ scales to zero, the Luttinger parameter Ku — > 1. 

The Hamiltonian (|16p is typically derived under the 
assumption that the electron-electron interactions are 
weak. On the other hand, it represents a stable low- 
energy fixed point of the theory, and thus should be valid 
beyond the weak-interaction approximation. Although 
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the nature of the low-energy fixed point can, in princi- 
ple, change at a finite value of the interaction strength, 
such a change would imply a quantum phase transition, 
which is generally not expected. The more likely sce- 
nario is that Hamiltonian (|16[) is the correct low-energy 
description of one-dimensional electron systems at arbi- 
trarily strong interactions. Under this assumption one 
should expect to be able to show that at low energies (i) 
the Hamiltonian ([9]) , ([4]) is equivalent to Eq. (fT6|) , and (ii) 
at the same time our expression (|14p for the electron de- 
struction operator transforms to Eq. pS]) . We now show 
that this is indeed the case. 



A. Low-energy Hamiltonian of strongly-interacting 
electrons in one dimension 

First we notice that the expressions ^ and (|16ap for 
Hp are very similar, as they both describe acoustic excita- 
tions in the charge channel propagating at speed Vp. Al- 
though it is natural to assume that K, 0, and 9 in Eq. ^ 
should be identified with Kp, (f>p, and 9p in Eq. (jf 6ap . re- 
spectively, this is not the case. The correct approach is 
to ensure that the physically observable quantities, such 
as the electron density ([3]) , have equivalent expressions in 
both theories. In the standard bosonization description 
based on Eq. (|15p the electron density is given by 



^]{x)^^ix)+i;l{x)ij^{x) = ^ + ^a,0p(x). (17) 

Comparing this expression with Eqs. pil)l . and pTjl . 
we conclude 



(t){x) = V2(f>p{x). 



(18a) 



Then, to preserve the proper commutation relations be- 
tween bosonic fields, one has to assume 



1 

3' 



(x) ^—9p{x). 



(18b) 



takes the form 

Ha 



H"" 



H'^y + H', 
I 



(20a) 
(20b) 

(20c) 



Thus the Heiscnberg model (|4]) is reduced to the tight- 
binding model of spinless fermions with repulsive inter- 
actions between particles at the nearest-neighbor sites. 

The steps leading from Eq. (g]) to (|^ are exact, and 
the spectra of the two Hamiltonians are identical at all 
energies. At energies much smaller than J one can sim- 
plify the Hamiltonian (|20p by bosonizing the Jordan- 
Wigner fermions. One starts by considering the non- 
interacting model given by Eq. (j20bp . The spectrum of 
that Hamiltonian is obtained as the sum of energies e{q) 
of independent spinless fermions. 



e{q) = J cos ql, 



(21) 



where the wavevector q varies from to 27r. 

In the absence of external magnetic field, one expects 
(5f ) = 0. According to Eq. ^ the band ^ is half- 
filled, and the two Fermi points are at 



77' QR = -TT- 



(22) 



The bosonization is accomplished by presenting the op- 
erator ai as a sum of operators destroying the right- and 
left-moving particles. 



ai = anil) + ql^I), 



where 



MO 



(23) 

(24a) 
(24b) 



Finally, substituting Eqs. (|18ap and (jlSbp into the Hamil- 
tonian ([9]), we recover Eq. (|16ap if 



K = 2Ko 



(18c) 



Turning to the spin Hamiltonian H„^ we note that 
there is a well-known procedure^^ of bosonization of the 
Heisenberg model One starts by converting spin op- 
erators Si to spinless fermion operators ai via the Jordan- 
Wigner transformation 

S'f = a\ai - i, Sf+iS'l = a\ exp ^-k ^ a]a}^ . (19) 

In terms of the spinless fermions the Hamiltonian H^, 



Since the bosonization description concentrates on the 
range of momenta close to the Fermi points, the discrete 
site number I is replaced here with the continuous co- 
ordinate ^. (Unlike the coordinate x of electrons, ^ is 
dimensionless.) The chiral bosonic fields Lpn and Lp^ sat- 
isfy the commutation relations 

bfl(e),%^fl(e')] = -2m5{i-a, (25a) 
[pdO.d^.vdO] = 2^*<5(e-n, (25b) 
[^L(0,%^fl(e')] = 0. (25c) 



Upon the bosonization (p4)) the Hamiltonian (|20bp 
takes the form 



(26) 



6 



The next step of the bosonization procedure is to convert 
to nonchiral bosonic fields 



^ 7:['PL + fR), If = 7:['PL - fR)- (27) 

As a result the Hamiltonian (l26l) takes the form 



Ho 



hv 

2^ 



(28a) 



with V = J/h and K, = 1. 

When the interaction term (|20cp is added to the Hamil- 
tonian (I20bp . and the bosonization transformation ((24|) 
is applied, the parameters v and IC change, and an addi- 
tional term appears in the Hamiltonian, 



V = 



25 



(27r5)2 



cos d^, 



(28b) 



with g ^ I. 

The Hamiltonian ()28p is equivalent to the spin part 
(jl6b|) of the Hamiltonian of weakly-interacting electron 
system, if one assumes 





= ■^'/'<T(C/n), 


(29a) 


m 




(29b) 


V 


= Van, 


(29c) 


IC 


Ka 
2 ■ 


(29d) 



In the absence of magnetic field, as the cosine term (j28bp 
scales to zero at low energies, the Luttinger parameter 
K, approaches 1/2, as required by the SU(2) symmetry 
of the problemji^ Thus Eq. (j29dp is consistent with the 
similar requirement 1 in the Hamiltonian (jl6bp . 



B. Bosonization of the operators Zi 
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In order to demonstrate that in the regime of low ener- 
gies |e| <C J the electron destruction operator (fT4|) takes 
the standard form one needs to bosonize the oper- 
ator Zi,^ in Eq. (fT4)) . We start with Zq^i. This operator 
acts on an arbitrary state in the Hilbert space of the 
Hamiltonian @ with N sites. If the spin at the site 
Z = is I, the operator Zo,x removes that site from the 
spin chain; if the spin is ti the outcome is zero. The 
Jordan- Wigner transformation defines one spinless 
fermion per each site with spin |. Thus the operator Zqj 
removes site I = from the tight-binding model (|20bp 
without changing the number of fermions Nj. 

To derive the bosonized form of Z^^i it is convenient 
to consider its effect upon the eigenstates of the non- 
interacting model (|20bp . The latter are Slater deter- 
minants of plane waves with wavevectors in the range 
< (7j < 27r and energies e{qj) given by Eq. (pij) . To 
determine the allowed values of wavevectors qj we as- 
sume periodic boundary conditions on the spin chain, 



(a) ZI 



2n q 



(b)ei 




FIG. 2: Ground state of the system of noninteracting fermions 
(|20bl) for the case of 16 sites and 8 fermions. Filled cir- 
cles show wavevectors and energies of the 8 occupied single- 
particle states, (b) Ground state of the same problem with 
15 sites. 



Sq — Sn- Because of the Jordan- Wigner string in the 
definition (jl9|) of the spinless fermions, their respective 
boundary conditions are either periodic, or antiperiodic, 
depending on the parity of their number Nf, 

«!v = i-^f'-'^^l (30) 
Thus the wavevectors of the fermions take the values 



1j 



N 



for odd Nf, 
f(j + i), for even iVj, 



(31) 



where j = 0, 1, . . . , — 1. 

The ground state of the Hamiltonian (j20bp is illus- 
trated in Fig. m^a). Upon removal of one site from 
the chain the allowed values of the wavevector increase 
slightly. 



N 



q, qj 



N- 1 



(32) 



The effect of the operator Zq.x on the eigenstates of 77^^ 
can be interpreted as follows. By removing a site from 
the spin chain, Zq^i creates a scattering potential for the 
Jordan- Wigner fermions near ^ = 0. Only forward scat- 
tering is present, and the wavcfunctions of the fermions 
acquire phase shifts proportional to the wavevectors. 



-if sgn (?) 



(33) 
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When the periodic boundary conditions are imposed on 
the fermions, the phase shift in Eq. ([55]) moves the al- 
lowed values of the wavevectors by q/N , in agreement 
with Eq. ([5^ . 

In the bosonization treatment of the spin chain one 
concentrates on the vicinities of the two Fermi points 
and g^, where the fermions are classified as cither 
right- or left-moving, Eq. According to Eq. 

the operator Zo,| transforms the fermion operators as 



1(0 



(34) 

Using the bosonized representation (I24p of the fermion 
operators and the commutation relations (|25p . we con- 
clude 



Zod = cxp [gL¥'L(0) + qRipni 



(35) 



Here we omit a numerical prefactor, which depends on 
the specific cut-off procedure used in the bosonization 
scheme, but can be considered to be of order unity. 

To find the bosonized expression for Zi^i away from the 
point Z = it is not sufficient to replace the arguments 
of (/?z, and ipR with I. Indeed, our derivation of Eq. ((35|) 
allowed for an arbitrary phase factor, which may depend 
on I. To determine this phase factor, we notice that 



iQl 



(36) 



where Q is the operator of the total momentum of the 
system. (This relation becomes clear if one notices that 
the operator e**^' shifts the spin chain by I sites to the 
left.) Summing the changes of wavevectors ([5^ for all 
particles between the Fermi points ((22)) . we conclude Zq^^t 
increases Q by 7r/2. Thus we obtain 



'-'1,1 = exp 



(37) 



where we have also substituted into Eq. ([35| the values 
of the Fermi wavevectors (j22l) . 

The apparent asymmetry between the left- and right- 
movers in Eq. (|37p can be understood by noticing that 
the interpretation of the changes of wavevectors of 
the fermions is not unique. Instead of assuming that as 
we remove a site from the spin chain, the fermion states 
in Fig. dja) transform into those of Fig. [D(b) by shifting 
to the right, Eq. (15^ . one can assume that all the states 
move to the left. 



N 



27r 



qj qj 



N ~1 N -1 



1j 



1j 



2tt 



N 



(38) 



As a result, the system arrives at a new ground state on 
the {N — l)-site lattice, which is the mirror image of the 
state shown in Fig.l^^b). One can repeat the above argu- 
ments leading to Eq. psp and obtain the new expression 



by replacing qR^L qR,L — 27r. Similar to the presenta- 
tion of the fermion operators a/ as a sum of two chiral 
contributions (l23ll. we conclude 



exp 



' exp 



. TT 3i i 



(39) 



Apart from these two contributions, the operator Zi^^ 
may contain terms corresponding to greater shifts of the 
fermion states in the momentum space. They can be 
obtained by adding any multiples of 2tt/N to the right- 
hand side of Eq. ([32]) . Such terms can be viewed as op- 



erator ([35]) combined with {ai}aR)"^ or (aij^ai)™. At 
low energies such terms are less relevant than the leading 
contributions and can be neglected. 

A similar bosonization procedure can be performed 
with operator Zi^^. In addition to the shifts ((5^ of 
the wavevectors of the allowed fermion states caused 
by the change of system size N N — 1, one also 
needs to account for the change in the fermion number, 
Nf — > Nf — 1. The latter changes the wavevector quanti- 
zation conditions (|31|) and removes a particle from either 
right or left Fermi point. The two contributions resulting 
from such treatment are 



exp 



' exp 



* 2 T^-^*- 



3i 



(40) 



Replacing the chiral bosonic fields in Eqs. (|39p and ([^01) 
with their non-chiral versions (l27l) , we find 



where the upper and lower signs correspond to 7 = t 
and J,, respectively. A generalization of the bosonization 
rule (|4T|) to the case of non-vanishing magnetization is 
discussed in Appendix IB] 



C. Two-step bosonization procedure for the 
electron operators 

In this paper we consider one-dimensional electron sys- 
tems with strong repulsive interactions, when the spin 
exchange between electrons is strongly suppressed, J ^ 
Ep. In such systems the bosonization of electron oper- 
ators ip-f{x) can be performed in two steps. At energies 
\e\ below the Fermi energy Ep the charge excitations can 
be bosonized, and the fermion operators take the form 
pi)) . This expression does not assume a specific rela- 
tion between e and J, so the spin excitations are ac- 
counted for accurately by the operators Zi^^. On the 
other hand, if |e| ^ J, the spin excitations can also be 
bosonized, Eq. (|4ip. To compare the resulting expression 



8 



for the electron destruction operators with those used in 
the standard bosonization procedure, Eq. ([TS]), we sub- 
stitute Eq. (jHJ) into Eq. (fTl]) . This substitution results in 
4 terms in the expression for ipj (x) . To identify the anni- 
hilation operator for the right-moving electron (jl5a[) , wc 
combine the first term in the brackets in Eq. (fT4|) with 
the first term in the right-hand side of Eq. (|¥T|l . This 
yields 



£ ^^lk'],x+<j>(x)] 



xg-«f;g±j[¥>(/)-ii?(0] 



g 2 F 



V27 



g-iS{x)^^4,{x)^±i[<p{nx)-^i>{nx)] 



(42) 



Expressing the bosonic fields via ^p.o- and 0p^a- with the 



help of Eqs. HH]) and ([^5]) . we find that our result 
is equivalent to the standard expression (jl5ap . Similarly, 
combining the second term in the brackets in Eq. 
with the second term in the right-hand side of Eq. (|4T|l . 
one reproduces the bosonized expression for the annihi- 
lation operator (|15bp of the left-moving electrons. 

To understand the meaning of the remaining two con- 
tributions to (x) , let us consider the momenta of the 
charge and spin excitations in Eq. By destroying 

the right-moving holon, we change the momentum of the 
system by kp = 2kp. Thus to obtain expression for an- 
nihilation operator of electron with momentum near the 
right Fermi point kp^ we chose the left- moving compo- 
nent of the operator Zi^^, which reduces the momentum 
change by fc^- By choosing the other component of Zi^^ 
we increase the total momentum change to 3k p. The 
physical meaning of such process amounts to removing 
an electron from the right Fermi point with simultane- 
ous transfer of another electron from the right to the 
left Fermi point. In interacting electron systems such 
processes are possible, but the resulting "shadow band" 
features tend to be weak. 



IV. GREEN'S FUNCTIONS 

A number of important physical properties of one- 
dimensional electron systems, such as the tunneling den- 
sity of states and the spectral functions, are expressed in 
terms of single-electron Green's functions, 

G+{x,t) - (ij^ix,t)iPl{0,0)), (43a) 
G^{x,t) = (V4(0,0)V^(x,t)). (43b) 

It is well known that in the limit |^| ^ oo the Green's 
functions show non-trivial power law behavior, with ex- 
ponents depending on the interaction strength. This 
behavior is easily obtained^ in the bosonization ap- 
proach based on Eq. In the case of strongly in- 
teracting electrons these results are valid at \t\ ^ h/J 



and adequately describe the physics of the system at low 
energies \e\ <^ J <^ Ep. 

In this paper we are primarily interested in the regime 
of intermediate energies, J ^ |e| ^ Ep. To find the 
Green's functions in this case, instead of Eq. ([TS]) one 
can use the more general expression ([H]). We start by 
transforming the electron annihilation operator (|14p to a 
more convenient form. 



Annihilation operator for strongly interacting 
electrons 



We first rewrite the lattice operators Zi^^ in Eq. ([T 
in terms of their Fourier components. 



J-TT 

oo 

.^(q) = J2 Zi,ye-^'^K 



(44) 
(45) 



Straightforward substitution of Eq. (j44|) into (fTi]) yields 



The two terms in the integrand correspond to removal of 
right- and left-moving holons, respectively. 

Expression (|46p can be simplified by noticing that 
Zj{q) is a 27r-periodic function of g. 



ip^{x) 



v27ra 



3tt 



^z^(g)e^(i+*)[*^--+*(-)l. (47) 



This presentation of the fermion operator is not entirely 
satisfactory, because the limits of integration here and 
in Eq. (|44p were set rather arbitrarily. Indeed, the argu- 
ment of the operator z^y (q) is the change of momentum of 
the spin subsystem when an electron is removed. Since 
spins in the Hamiltonian @i are attached to lattice sites, 
the momentum changes by g or g ± 27r, q ± An, etc. are 
equivalent. This symmetry is lost in the bosonized ex- 
pression ([Tf|) , but can be restored by extending the limits 
of g-integration. 



-i0(x) 



V2^ 



dq 
2^ 



)[fcJ-2;+0(x)] 



(48) 



The origin of the ambiguity in the definition of the 
fermion operators can be traced back to the bosonization 
of the holon operators in Eq. ((8b)) . Indeed, by definition 
([7]), the operator l{x) has only integer eigenvalues equal 
to the number of electrons in the region of space from —00 
to x, whereas its bosonized expression l{x) — —^\x -V 
0(a;)] does not explicitly possess this property. To enforce 
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the discreteness of charge in Eq. IM]) one can understand 



as 



;=i:[fej,2;+0(a;)] 



Z,.^(5(i[fc^x + 0(:r)]-/). (49) 



bosonic fieldsi^ In the most important regime |t| ^ U/T 
one finds 



X — Vpt ±ia J \x + Vpt =F za 



(53) 



Then, upon the Fourier transformation (|44p , one recovers 
Eq. (gSl). 

A similar procedure of bosonization of fcrmions while 
preserving the discreteness of their number was suggested 
by Haldane.'^^ Apart from the two terms corresponding 
to the right- and left-moving fermions, the expression 
for the fermion operator contains multi-particle contribu- 
tions with wavevectors near ±3fci?, ±5kp, etc. In a typ- 
ical bosonization calculation these additional terms give 
much smaller contributions than the leading ones. Thus 
the difference between the results obtained using the two 
bosonization schemes is smaller than the accuracy of the 
bosonization approximation, and can be ignored. In our 
case, by using expression (|47p instead of ([35]), one obtains 
small spurious features in the spectral function originat- 
ing from the arbitrarily chosen integration limits. Thus 
from now on we use Eq. 



where 



B. Green's functions at intermediate energies 

At small J the time evolution of the spin degrees of 
freedom is very slow, with the typical time scales of order 
h/J. Therefore to find the Green's functions describing 
physical phenomena at relatively high energies \e\ 3> J, 
one can neglect the time-dependence of the correlators 
of operators Zi^-y. Then the substitution of Eq. ([^5]) into 
the definitions ([43]) of the Green's functions yields 



G±(x,t) 



g^{x,t). (50) 



Here the static correlators (q) are defined by 

c+(q) = ^(Z,,4^)e-'^'', (51a) 
I 

c-{q) = ^(4^Z,,)e-'^'', (51b) 



where (...) denotes thermal averaging over the equilib- 
rium states of the Hamiltonian ^ . The time dependence 
is contained in the correlators g^, defined in terms of the 
charge variables, 

g+(x,t) = (^e^li^+i)'Pi^^t)-H^,t)] (.-m+i)'Pio,o)-o{o,o)]'^^ 
' (52a) 

(52b) 

These correlation functions are easily computed using 
the standard techniques for averaging the exponentials of 



(54) 



The correlators (g) are determined by the properties 
of the Heisenberg spin chain Unlike the holon cor- 
relators g^{x,t), in general c^{q) cannot be computed 
using the bosonization approach. We discuss their be- 
havior in detail below. 



C. Static correlators c^{q) 

To find c^{q) one has to perform the averaging in 
Eqs. ([5T|) over the eigenstates of the Heisenberg Hamil- 
tonian The results depend crucially on the relation 
between the temperature T and the exchange constant 
J. At T ~ J it can only be studied numerically. Some 
analytical results can be found in the cases of high and 
low temperatures. 



1. General mathematical properties of (q) and (q) 

We start by establishing interesting relations between 
the functions c^{q) and c~{q), which follow from their 
definition (j5ip and do not depend on the temperature 
or the specific form ^ of the Hamiltonian of the spin 
chain. We first notice that the correlators (Z; ^Z^^) and 
(Zn ^,Z, „,) are real, and satisfy the following relations 



(■2'o,7-^0,7) 
(■^i,7-^0,7) 
(^0,7'^'-7) 



= 1, 



= LZ_ 



(^0,7^-i,7) 



From the definitions ([5T|) it then follows that the corre- 
lators c^{q) are real and even functions of q. Also, the 
definition (j51ap of (q) can be rewritten as 



l + 2Re£(Z;,^4^)e- 



■iql 



(55) 



1=1 



Furthermore, one can establish a simple relation between 
operators Zi^Z^^ and Zq^Z^^. Both of them add and 
remove a site with spin 7 at different positions on the 
spin chain. If a site is added first, the numbering of all 
the subsequent sites is shifted by 1. Thus Zj ^Zq — 
Zq ^Zi_^ ^. (Commutation of operators Z^ ^ is discussed 
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in more detail in Appendix [X]) Then, using Eq. (|51bD 
one obtains 

Substituting this relation into Eq. (|55p one finds 

/•^ fin' sin 

c+{q) = l + c-{q)cosq + J ^^—^c-{q'). (56) 

Similarly, one can express (q) in terms of c+ (g) , 

-/X -I-/ \ r dq' q — q', , 

S (9) = S ((7)cos(7+ / — cot — ^ smg c:J:((7 ). 

(57) 

The integrals over q' in Eqs. ([55]) and ([Ff|) should be 
understood as the principal value. 

The relations and ([57)1 simplify considerably at 
q = and q — tt. In these cases wc find 

c+(0)-c;(0) - l + (4^Zo,^), (58a) 
c+(7r)+c-(7r) = l-(4^Zo,^). (58b) 

The operator Zq ^Zq ^ destroys and then recreates a site 

with a given spin 7. Thus the average (Zq^Zq^) is the 
probability to find spin 7 at site / = 0. In the absence of 
magnetic field one expects {Zq .^Zq .^) — i, resulting in 

c+(0)-c;(0) = ^, (59a) 

c+{n)+c-{7T) = 1. (59b) 

The relations and (|57p follow from the definition 
of the correlators (q) , and are not sensitive to the spe- 
cific Hamiltonian ([4]) of the spin chain. The only prop- 
erty of the Hamiltonian ^ important for the relations 
([55)) is the spin symmetry. At the same level of uni- 
versality one can find the full dependences c^{q) in the 
high-temperature limit. 

2. High temperature T 'S> J 

At high temperatures, correlators similar to c^{q) have 
been studied by Penc and Serhan.'^^ They noticed that at 
J /T — > the spins are completely uncorrelated, and the 
correlator (Z^ ^Z^ ^) is simply the probability of finding 
|/| random spins at sites 0, 1, . . . Z— 1 (or /, / + 1, . . . , —1 for 
negative I) in a given state 7. Thus (Z; ^^9,7) — 1/2'' 
and, similarly, {Zq .^Zi^) — l/2l'l+-^. Substituting these 
expressions into Eq. ()5ip . one finds 

c+{q) = 2c-{q)^—^ , T » J. (60) 

' ' 5 — 4 cos q 

One can check explicitly that these results are consistent 
with the relations (|56p and ([57)1 . 



3. Zero temperature 

At zero temperature the correlators c^{q) and c^{q) 
carry non-trivial information about the spin correlations 
in the ground state of the antiferromagnetic Heisenberg 
spin chain ([4]). Although the model ([4]) is exactly solv- 
able, no exact results are known for the correlators (q) . 
The quantity analogous to c~{q) was first studied by 
Sorella and Parol^S, who used the results of numeri- 
cal diagonalization of spin chains of up to 22 sites4i 
Their results indicated that c~{q) is extremely small at 
< (7 < 7r/2, whereas C^lq) ~ 1. They interpreted this 
behavior as the effect of "spinon pseudo Fermi surface," 
which can be rephrased as follows. 

Let us consider the ground state of the fermionized 
version ([^0)1 of the Heisenberg model. Ignoring the inter- 
actions ()20cp between the fermions, one can picture the 
ground state as shown in Fig.[2ja). The function [q) is 

defined in terms of the correlators (Zq ^Zi ^). By remov- 
ing a spin-f site, the operator destroys a fermion. 
Ignoring for the moment other aspects of one con- 
cludes that c^{q) should vanish at — 7r/2 < q < tt/2, as 
the fermion states with those values of the wavevector 
are empty, see Fig.[2](a). A similar argument for the cor- 
relator c^{q) shows that it should vanish in the region 
of g-space below the Fermi surface. Noticing that in the 
absence of magnetic field specific spin direction is unim- 
portant, one concludes 

c+{q) = at|<g<^ (61a) 
c-{q) = at-|<g<| (61b) 

Conditions ([CT)) refer to the values of q between — 7r/2 and 
37r/2; outside of that region they can be inferred using 
the 27r-periodicity of the functions c^{q)- 

The above picture neglects two important aspects of 
the problem: (i) the fermions interact with each other, 
Eq. ()20cp . and (ii) the operator Z;^ not only destroys 
a fermion, but also removes a site from the spin chain. 
There is no a priori reason to expect that the effects 
of these approximations should be small. However, the 
bosonization approach to calculation of c^{q) near the 
Fermi points q — ±tt/2 (outlined below) shows that al- 
though each effect is significant, they mostly compensate 
each other. Numerically, this compensation result a^^i^^ 
in c;^((7) < 0.01 for \q\ < 7r/2. Based on their numeri- 
cal results, Sorella and Parola'*° conjectured that in an 
infinitely long spin chain c^{q) = at |q| < 7r/2. From 
this point of view, the numerically small values of c~{q) 
in this interval should be viewed as a finite-size effect. 

A more detailed study of the static correlators of the 
spin chain ((D) was performed by Penc et al.^^ They stud- 
ied both correlators cf^{q) and c~{q) (in slightly differ- 
ent notations) and found that not only c^{q) is numer- 
ically small at0<g<7r/2, but also c^(g) is small at 
7r/2 < g < TT. Although these observations are consistent 
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with the single-particle Fermi surface prediction ((6T|) , no 
significant size effect was observed. Penc et al?^ con- 
cluded that contrary to Eq. ((6T|) . correlators c+((7) and 
c~ (q) do not vanish exactly in the "forbidden" regions of 
the wavevector q. 

We now show analytically that although the conditions 
(j6ip may be a good approximation of actual behavior 
of correlators {q) , they contradict to the spin-rotation 
symmetry of the problem. We showed in Sec lIVCTl that 
the correlators c^{q) and c^{q) are related to each other, 
and if one of them is known on the interval of q of length 
2tt, the other can be obtained using the relations ([55]) 
and ([ST]) . Conditions ((6T|) define c+ (q) over half the pe- 
riod 2tt and c~{q) over the other half. This is sufficient 
to uniquely define both functions at the remaining re- 
gions of (?-space. To accomplish that one can solve the 
integral equation ([57)1 with conditions ([CTjl . This mathe- 
matical problem belongs to the class of singular integral 
equations, and can be solved by applying the well-known 
techniques;^ see Appendix O The solution has the form 




a Q 9 Q o — I — S g B B B n b 

Jl/2 Jl 

Wavevector q 

FIG. 3: The solution of the integral equation (I56|l with con- 
ditions (|6Hl . The solid and dashed lines show the behavior of 
c^{q) and c~{q) described by Eqs. (|62a|l and (|62b[l . respec- 
tively. Squares and circles represent the results of numerical 
solution2flii^ for the 26-site chain Q. 



<(?) = 



Co 



exp 



^cosg 
for — 7r/2 < q < 7r/2 and 



Tr/2 



In 



Co 



-7r/2 



tt/2 



dq' 



(62a) 



\/\cosq\ 



exp 



Tr/2 



In sm — - — aq 



(62b) 

for 7r/2 < q < 37r/2. The normalization constant cq is 
found using Eq. ([5? 



1 



Co 



2 cosh ^' 



(63) 



where G « 0.91597 is the Catalan's constant. The solu- 
tion ((62a| . ((62b|) is plotted in Fig. H 

As expected, our solution gives an excellent approx- 
imation to the numerical data of Ref. [l^. Neverthe- 
less, the results ([6T|) . (|62)) are not exact. Indeed, from 
Eqs. (Pa|l and ((62b)l we find 



c+(0) 
c~{t:) 



1 + e-iG/TT 
2 



1.525, 
0.475. 



(64a) 
(64b) 

= 0, 



Taking into consideration that c+(7r) — c;^(0) 
Eq. (pT|) . we find that the spin-symmetry conditions ^ 
are not satisfied exactly. 

To gain better insight into the properties of the corre- 
lators c^{q), let us consider the case when the wavevector 
approaches the Fermi point, g — > 7r/2. In this regime one 
can apply the bosonization approach of Sec. IIII Bl Using 
the expressions (j4ip for the operators Zi^^ and applying 
standard techniques'^ for averaging the exponentials of 



bosonic field over the ground state of the Hamiltonian 
(I28a|) . we find 



(^(,7^0,7) 



a + il 



1/2 



a ~ il 



52 xlCc+^-i) 



a2 + ;2 



1/2- 



(65) 



In the limit of low energies the Luttinger-liquid constant 
JC 1/2. Using this value one can find the correlator 
c^{q) near the Fermi point q = n/2 by substituting the 
component of the Green's function (|65p given by the first 
term in the brackets into Eq. (j51ap and replacing the sum 
over I with an integral: 



C+(<Z) = 



a + il 



1/2 



-^(q-^/2)l^l_ (66) 



It is important to note that the integrand is analytic 
in the lower complex half-plane. Thus cj^{q) = for 
q > tt/2. This conclusion agrees with the prediction of 
the spinon Fermi surface approximation (j6ip . At g < 7r/2 
the integration in Eq. (|66p is straightforward, and we find 



ef- 



Cj{q) = X- 



(67a) 



where Q{x) is the unit step function, and the value of 
the numerical coefficient % ~ 1 cannot be determined 
within the bosonization approach. A similar calculation 



for [q) results in 



c^ (<?) = X 



(67b) 
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also in agreement with the prediction ([61]) of the spinon 
Fermi surface picture. 

The inverse-square-root singularities ([67]) of the corre- 
lators c^{q) at the Fermi point are consistent with the 



numerical data of Ref. l20| and our expressions (|62p . see 
also Fig. [3l In particular, the asymptotes of the expres- 
sions (p^ at q ~ n/2 are given by Eq. ([57)1 with x ~ 0.85. 

It is important to note that the results ([S7)) of the 
bosonization treatment of the static correlators c^{q) do 
not prove the validity of the spinon Fermi surface picture. 
Indeed, the bosonized expressions (|4ip for operators ^ 
are only applicable asymptotically near the Fermi point, 
q —^ ±7r/2. Consequently, the presence of the step func- 
tion 8 in the results ([57)1 should be interpreted as 



lim 



(68) 



Let us now show that the correlator c^{q) not only does 
not vanish at g > n/2, but in fact diverges at<7^7r/2+0. 

We start by noticing that the conclusion c:|(q) = at 
q > 7r/2 holds only for /C = 1/2. Indeed, at /C ^ 1/2 
the last factor in the expression (|65p is no longer unity. 
More importantly, it is no longer analytic in the lower 
complex half-plane, leading to the non- vanishing value of 
the integral (|66p at q > tt/2. 

The low-energy properties of the Heisenberg spin chain 
(III) are adequately described by the sine-Gordon Hamil- 
tonian ()28|) . The parameters /C and g renormalize at low 
energies or, equivalently, long length scales L as 



(69) 
(70) 
(71) 



Thus at low energies K. deviates slightly from 1/2. This 
deviation gives a correction to the correlator (j66p which 
can be obtained by expansion of Eq. (|65p in powers of 
y = 4(/C-l/2), 



JC = 


1 y 

2 4' 


9 = 


TTVy, 

1 


y = 





2 greater than Eq. ((73|) . [Apart from direct calculation, 
this can be shown to follow from the spin-rotation sym- 
metry of the problem, see Appendix |DJ] We therefore 
conclude that at q slightly above the Fermi point tt/2 
the correlator c+ is given by 



32 



TT 



(74) 



The above calculation was performed to second-order 
perturbation theory in coupling constant y and did not 
account for its scaling ([7T|) . Since the scaling is logarith- 
mic, and thus slow compared to the leading power-law 
behavior in Eq. ()74p . one can simply substitute the ex- 
pression ([7T|) into (fZl)) choosing the proper value of the 
length scale L ~ l/{q — tt/2). We therefore conclude 



that near the Fermi point q 
behaves as 



Tr/2 the correlator c^{q) 



<(g) 



3tt 1 

32 ln^(g- 



q 



- - 

2 



0. 



(75a) 



Analogous calculation for the correlator gives 
-7=, 1 



f + 0, 



S (5') = 



32 ln2(f-5) 



0. 



(75b) 



Our results (l75l) reaffirm our earlier observation that 



the correlators c^{q) do not vanish exactly on one side of 
the Fermi point, as expected from the spinon Fermi sur- 
face picture, cf. Eq. ([61]) . Instead the correlators diverge 
at q tt/2, albeit slower than on the "main" side of 
the Fermi point. The numerical data of Ref. [10 shown in 
Fig. [3] indicate that the small values of €+((?) at g > 7r/2 
do increase near the Fermi point. However, studies of 
much longer systems are needed to verify the asymptotes 



a 



a + il 



1/2 



In 



(72) 

At q < tt/2 the correction ([7^ is small compared to the 
leading term (j67ap and can be ignored. However, unlike 
Eq. (|67ap . the correction ([72]) does not vanish at q > tt/2, 



Sc+iq) 



32 



(73) 



where x is the same numerical coefficient as in Eq. ([57]). 

The deviation of the Luttinger liquid parameter /C from 
the limiting value of 1/2 is not the only source of correc- 
tions to c^{q)- Additional corrections originate from the 
small sine-Gordon term (|28bp . One can account for this 
term to second order perturbation theory using the stan- 
dard techniques. The resulting correction is factor of 



SPECTRAL FUNCTIONS 



The spectral functions {k, uj) and {k, uj) of a one- 
dimensional electron system are defined as Fourier trans- 
forms 



A^{k,L0) 



-ikx-\-iiut / 



G^ix,t) (76) 



of the Green's functions ([43]) . The components A+(fc,w) 
and A~(k,Lu) characterize the particle and hole parts of 
the excitation spectrum, respectively. In this section we 
study the behavior of the spectral functions at frequen- 
cies UJ in the range 



(77) 
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Under these restrictions one can express the Green's func- 
tions in the form ([50]) . ([STjl . ([55]) . Performing the inte- 
gration with respect to x and t, one finds 





[q-qk + qujV~'^^ [qk + qu>- qV^'^" 



where 



(78) 



(79) 



According to Eq. (|75|) at small frequencies u; <C 
the spinon wavevector q must be close to qk. This can 
be understood in terms of the energy and momentum 
conservation laws. Since the energy u is small, the holon 
component of the electron must have momentum near 
±kp. If the electron momentum k is near one of those 
values, the spinous carry no momentum, i.e., q = (or 
27r), in agreement with g sa « 0. On the other hand, if 
k is not near ±kp, the difference of the momenta k—kp — 
nqk is transferred into the spin subsystem. (Note that the 
electron density n = kp/n.) 



A. Zero momentum peak 



This result reproduces the expressions for the tunneling 
density of states at e/ Dp reported in Ref. [Slj. 

It is important to note that the Gaussian peak 
in the spectral function (|80p cannot be obtained 
within the standard Luttinger-liquid theory of spectral 
functions, '^'^'^^ which applies only near k — ztkp. The 
peak ((80)) is due to the holon states with wavevectors 
near kp — 2kp, and the spin excitations with q near tt, 
well below the spinon Fermi surface, see Fig. O 

An interesting consequence of the fact that the Gaus- 
sian peak is dominated by spin excitations away from 
the Fermi points is the dramatic difference of the peak 
heights for the spectral functions and A~ . In- 
deed, as we discussed in Sec. IIVCI at zero temperature 
c+(7r) <C c~{tt), whereas at T » J one has c+ = 2c~ . 
Thus the heights of the peaks at fc = depend strongly 
on temperature, and the peak in A~ is much more pro- 
nounced than that in at T ^ J. 

In our derivation of the expression (|5D|) for the spec- 
tral functions we replaced q — > in the arguments of the 
spin correlators (q) in Eq. ([75]) . This procedure is well 
justified if qk is not too close to points 7r(2s — l)/2 (with 
s = 0,±l,±2,...), where at zero temperature the corre- 
lators c^{q) have sharp inverse-square-root singularities, 
see Fig. [3] In addition, one can still use Eq. ([50)1 if the 
singularities of c^{q) are smeared by finite temperature 
by Sq ^ T/J ^ qi^. However, at very low temperatures 



To find the momentum dependence of the spectral 
functions (j78p at low frequencies, one can replace q —^ qk 
in c^{q) and The remaining integration is straight- 
forward and gives 



(fc,w) 



e(±^)c±(gfe) / D 



2«^4r(c(fc)) \h\L^ 



where 



i-C(fc) 



(80) 



T< 



TlLjJ 

"d7 



(83) 



the approximation leading to Eq. (j80p fails near qk = 
n{2s — l)/2. Thus at low temperatures the spectral 
functions have non-trivial behavior in the vicinity of 
k = (2s + l)kF- Below we consider the zero-temperature 
behavior of the spectral functions near these points. 



c(fc) = c++c = 



K 



Uh 

h,p 



1 

2K' 



(81) 

In the low-frequency limit h\uj\/Dp — s- the spectral 
function A^ has a sharp Gaussian peak as a function 
of the electron momentum centered at fc = 0. In the 
spin-incoherent limit J /T — > a similar Gaussian peak 
in the momentum dependence of the spectral function 
was found by Fiete et al.'^^ Our expression ([80]) is valid 
at arbitrary J/T, with the temperature dependence en- 
tering Eq. (I80p via the functions c^{q). In particular, 
using the numerical results2£ for (q) we can access the 
zero-temprature limit. 

The peak (|80p gives the leading contribution to the 
density of states at low energies 



dk 
2TTn 



A^{k,e/n) 




In 



(82) 



B. Fermi surface features 

When the electron wavevector approaches the Fermi 
point kp we have qk — 7r/2. From Eqs. (j67ap and 
the symmetry c^{q) = c^{—q) we conclude that near 
q — — 7r/2 one can approximate c+(g) as 



(84) 



kl 



To explore the fine structure of the Fermi surface features 
at zero temperature, instead of the approximate expres- 
sion ([50)1 one should use the more general formula (|78p . 
Then instead of one singularity at qk — —tt/2 one finds 
two singularities at qk = — 7r/2 ± q,^. Given the defini- 
tions ([75)) . one concludes that the spectral function 
has singularities at two points k = kp±uj/vp, above and 
below kp. 

In evaluating the singular terms in the spectral func- 
tion (I78p one can approximate the exponents C,^ by their 
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values at q = —Tr/2. Using the expressions (|54|) we then 
find = C-yr/2 ^ = Co, where 



Co 



{2-Kf (1 



16K 



8K„ 



(85) 



The stronger of the two singularities is the one above 
kp, at k = kp+u/vp. It appears when the inverse square 
root singularity (|84p in c^{q) is near the lower limit of 
the integral ([75)1 . The most singular term in the spectral 
function can then be found by extending the upper 
limit to infinity, resulting in 



A+ik,L0) 



27rr(Co + i)r(Co) V 2"pfcF 



X [2a;|w — Vp{k — k[ 



(86) 



where the constant B is defined in terms of the beta 
function 

r B(i-Co,5), uj-Vpik-kF)^+0, 
B = { (87) 
iB(i-Co,Co), Lo^Vp{k~kF)^ -Q. 



The power-law singularity (|86p of the spectral function 
A^^ at fc = kp+Lu/vp with the exponent Cq — 1/2 is consis- 
tent with the results of the Luttinger-liquid theory4^i^ 
It is worth noting that at 



K 



2Kp>6 



4\/2« 0.343 



(88) 



the spectral function ((86)) diverges at k — kp + Lu/vp. 
Nevertheless, the dominant contribution to the density 
of states (|82l) at low energies is given by the Gaussian 
peak dHO]) at fc = 0. 

The second singularity of the spectral function is 
below kp, at k = kp — uj/vp. It emerges when the sin- 
gularity ([M)) of c+(g) is at the upper limit of the inte- 
gral ([78]) . As k approaches kp — oj/vp from above, the 
width of the integration regions shrinks to zero, and at 
k < kp ~ Lo/vp the integral vanishes. The singularity at 
k ^ kp — uj/vp has the form 



At{k,u;) 



9(fc - ^kp - oj/vp))xUc^l2vpf^°-''^ 

27r[r(Co + i)]'v/2^ 
y,{2Ljf°-^[uj + Vp{k~kp)]^''. (89) 



Unlike the feature (l86|) above kp, this term always van- 
ishes at the singularity because Co > 0. The power-law 
feature in the spectral function with the exponent Co was 
obtained earlier— in the framework of the Luttinger 
liquid theory. 

The fact that the power-law feature ((89|) appears only 
on one side of the point k — kp ~ uj jvp is a. consequence 
of the "spinon Fermi surface" approximation. In a more 
careful treatment, the correlator c^(q) does not vanish 
at g < — 7r/2, but instead has an inverse-square-root 
singularity with an additional factor l/ln^(g -I- 7r/2), cf. 



Eq. (|75ap . Thus the cusp ([89)1 should appear on both 
sides of the point fc = kp—uj/vp, albeit with an additional 
logarithmic suppression factor at fc < fci^ — u)/vp. In the 
Luttinger liquid theory this feature would be caused by 
the marginally irrelevant perturbation cos(2V2 0cr) that 
was not accounted for in Refs. |45||46| . 

At cij < the spectral function A~ has similar power- 
law singularities at fc = kp ± Lu/vp. The respective ex- 
pressions can be obtained by replacing uj ~lu and 
k — kp ^ kp — k m Eqs. (|86p and 



C. Shadow band features 

The singularities in the spectral functions emerge when 
excitations in both charge and spin subsystems are near 
the Fermi points kp and (2s — l)7r/2, respectively. In 
particular, if the holon momentum is kp = 2kp and the 
spinon momentum nq = (2fci?/7r)(7r/2) = kp, the elec- 
tron momentum is 3fci?. Formally, the features in the 
spectral functions are caused by the singularities (p7)) ap- 
pearing inside the narrow integration region in Eq. (|78p . 

Similarly to the singularities near kp, one expects to 
find two features, at fc = 3fci;' ± to/vp, when either of 
the limits of integration crosses the point q = Tr/2. Since 
the features in the spectral functions near 3k p are weaker 
than the ones near kp, we discuss only the stronger of the 
two singularities in A^ . It appears when the lower limit 
Qk ~ luj is near 7r/2 and corresponds to fc = Sfc^ + ^/vp. 

The behavior of the spectral function is controlled by 
the exponents ([5^ at q = 7r/2, which can be expressed 
in terms of a single parameter 



Ci 



{iK - 2) 
16^ 



{iKp 



8K, 



(90) 



Cg = Cq'^ ^ 3/2 = Ci- Then the singular term in the 
spectral function takes the form 



9(3fcF + uj/vp ~ k)x{a/2vpf^^+^'^ 

27r(Ci + i)[r(Ci)]2y2^ 
X (2ujf^+^/''[uj - vp{k - 3kp)]'^'-^^^. (91) 



AUk,Lo) 



Similarly to the singularity ([89)) . one will find a weaker 
feature on the other side of the singularity, i.e., at fc ^ 
(3fc_F +uj/vp) + 0, if instead of Eq. (|67ap a more accurate 
approximation (|75ap for the correlator c+((?) is applied. 

The features in the spectral functions at fc = 3fci? -I- 
uj/vp have been observed in numerical data for the 
infinite-?/ Hubbard model by Penc et al.,^^ who iden- 
tified it with the so-called shadow band.~ Our formula 
(|9ip provides analytic expression for the spectral func- 
tion at the shadow band position in the limit of low fre- 
quencies Lo <C Dp/h. Unlike the numerical treatment of 
Ref. I21I, our result is not limited to the Hubbard model 
with only on-site repulsion and can be applied to systems 
with any interaction range. 
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In addition to the shadow band feature near k = Sfcj?, 
the periodicity of the correlators c^{q) results in sin- 
gularities of the spectral functions at all odd multiples 
of kp. Similar to the features near kp and 3k p, one 
finds a pair of singularities at fc = (2s + l)kp ±Lu/vp for 
s — 2,3, . . .. The stronger singularity in each pair is the 
one at fc = (2s + l)kp + Lu/vp, where one finds 

A^{k,uj) oc 

x{tj-wp[fc- (2s + l)fcF]}'^'"^/^ (92) 

with 

_ [(2s + l)j^-2p _ [(2s + l)j^,-lF 
16K " 8Kp ■ 

At s = 0, 1 this expression is consistent with our earlier 
results ((86)) and ((9T|) . For a given s the strongest (inverse 
square-root) divergence of the spectral functions at fc = 
(2s -|- l)fcF + ^/vp is achieved when (s = 0, i.e., Kp = 
l/(2s+l). Since in the case of strong repulsion < 1/2, 
this condition cannot be satisfied for s = 0; the lowest 
possible value of (o is 1/16. 



VI. SUMMARY AND DISCUSSION OF THE 
RESULTS 

In this paper we have developed the theory of one- 
dimensional electron systems in the regime of very strong 
interactions. This regime emerges when the repulsion 
between electrons strongly suppresses exchange of their 
spins, J <ti Ep. Our theory is based upon the Hamil- 
tonian H = Hp + H^, with the charge part ^ brought 
to the form ^ by means of conventional bosonization 
while the spin contribution H^, is the Hamiltonian of the 
Heisenberg spin chain The most important ingre- 
dient of the theory is the expression ([5]) for the electron 
creation and annihilation operators in terms of the charge 
and spin degrees of freedom. 

In our technique the charge excitations are bosonized, 
and thus the applicability of the results is limited to en- 
ergies well below Ep. On the other hand, the spin exci- 
tations are treated more carefully, so we can access the 
energy scales both below and above J. At energies be- 
low J the standard approach based on the bosonization 
procedure and the Hamiltonian can be applied. 
We showed in Sec. |lll] that at e ^ J our expressions ([5]) 
for the electron operators reproduce the bosonization for- 
mulas ([15]). The advantage of our method is that unlike 
the bosonization procedure (jl5p it can also be applied at 
energy scales e ^ J. 

The main difficulty in applying our technique is the 
need to find the Green's functions of operators Zi^^ with 
the Heisenberg Hamiltonian ([4]). The problem is simpli- 
fied for the most interesting case e ^ J, when the slow 
time dependence of the spin degrees of freedom can be ig- 
nored. In this case the single-particle Green's functions of 



electrons can be expressed in terms of the static spin cor- 
relators c^{q), which were studied in Refs. [20ll23l39ll40l . 
Additional useful properties of these correlators are de- 
rived in Sec. HVCl 

We have applied our technique to the calculation of the 
spectral functions of strongly interacting one-dimensional 
electron systems in Sec. |Vl Experimentally, the spectral 
functions can in principle be studied by angle-resolved 
photoemission spectroscopy. However, we are not aware 
of such experiments on one-dimensional conductors in 
the regime of strong interactions. A more promising ap- 
proach is to observe momentum-resolved tunneling be- 
tween two parallel quantum wires. Such measurements 
have been recently reported by Auslaender et al^i^ 
The regime of strong interactions can be achieved in a 
quantum wire by reducing the electron density. Unfor- 
tunately, at low densities the effects of disorder are also 
amplified, and the tunneling into states of given momen- 
tum is no longer possible.'*^ Thus to observe the predicted 
features in the spectral functions one would have to man- 
ufacture quantum wires with even less disorder than in 
Refs. 48.49. 

Our results can be compared with previous theoretical 
studies of the spectral functions of one-dimensional elec- 
tron systems. In particular, Meden and Schonhammer— 
and Voit^ studied the spectral functions in the frame- 
work of the bosonization approach. In the case of a 
strongly interacting system their results are valid only 
at energies below J. However, they can still be com- 
pared with our results at w ^ J /Ti by taking the formal 
limit in the bosonization results, as s uggest ed by 

Pcnc et al?'^ The spectral functions of Refs. [45ll46l show 
singularities at k — kp ± uj/vp, which are not sensitive 
to the spin velocity Va- As expected, the power-law be- 
havior of the spectral functions at those singularities is 
in agreement with our results (j86p and (j89p . Similarly, 
in the Uo- ^ limit, the singularities at w = ±Va-{k — kp) 
of the spectral functions of Refs. |45||4^ show the same 
power-law frequency dependence as the Gaussian peak 
([5U)) in the tail region (at fc = fc^). 

Since the bosonization technique accounts only for the 
electrons near the Fermi points, our results for fc away 
from ±kp cannot be compared with those of Refs. I45ll46l . 
In particular, the main contribution to the tunneling den- 
sity of states at low energies is due to the peak in the 
spectral functions centered at fc = 0, Eq. (|5D|) . Conse- 
quently, our expression ([5^ for the density of states in 
the case of the Hubbard model {K = 1) gives a larger 
result v cx e~^/^ than the Fermi-surface contribution 
v (X e"'^/^ by Penc et alM- The inverse square-root de- 
pendence of the density of states on energy was obtained 
earlie r^^i^"^i^^ in the case of T ^ J. Our results show^ 
that the same dependence holds also for T < J. Phys- 
ically the enhancement of the density of states at low 
energy is analogous to that in the X-ray absorption edge 
problem,— with the spin excitations creating the the ef- 
fective core-hole potential for the holonsj ^^'^^ 

A detailed comparison can be made between our re- 
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suits for the spectral function and those of Penc et al^ 
The latter work studied numerically the quarter-filled 
Hubbard model in the limit of infinite on-site repulsion 
U. To compare their results with ours one should assume 
K = 1, and consequently Co = Ci = 1/16. Our results 
([86|) and (j9ip indicate that power-law peak with expo- 
nents Co,i ~ 1/2 = —7/16 should appear at k = kp+uj/vp 
and k — 3kp + uj/vp. The data of Ref. [Uj does show sin- 
gularities at those lines in the (fc,a;) plane. In addition, 
as w — *■ the singularity at k — kp + uj/vp is expected 
to grow as oj^^/^^, whereas the one at fc = 3kp +uj/vp is 
expected to be suppressed as upl^^ . The data of Ref. [U 
does show this qualitative behavior. Finally, the data^i 
clearly shows a weak feature at k — kp ~ oj/vp, which 
becomes more prominent at a; — )■ 0. This feature is con- 
sistent with our result (|89p which at Co = 1/16 behaves 

as W-15/16[fc-(fc;.-w/z;p)]l/16. 

One should note that the numerical data^^ does not 
show a peak at fc = that we expect based on Eq. ([80)1 . 
At Lo +0 the spectral function appears to be very 
small. This can be understood as a result of smallness 
of c+(7r) w 0.044. At w ^ -0 the spectral function^^ 
is not small, but instead of a Gaussian peak it shows a 
rather flat minimum at fc = 0. This can be understood by 
noticing that the fc-dependent prefactor c~((7A;)/r(C(^)) 
in Eq. (j80p has a minimum at fc = 0. In the limit 
\u)\/Dp the last factor in Eq. ((80)l dominates, and 
we find a peak. However, in a finite system one cannot 
access the values of ?i|a;| below the level spacing. Sub- 
stituting the parameters of the Hubbard chain used in 
Ref. Hi] into Eq. (|80| and using the approximation (|62bD 
for c~{q), we find 

\nA-(k,u}) = const + f 2.03 - 4- In | fc^ (94) 

' \ TT^ "-1^1/ 

at fc 0. Thus the spectral function should have a peak 
if \n{Dp/h\uj\) > 10. On the other hand the finite level 
spacing on the quarter-filled lattice of 228 sites^i limits 
the frequencies such that \n{Dp/h\cu\) < ln(228/27r) « 
3.6. Thus to find the peak at fc = significantly longer 
systems should be studied. In the spin-incoherent regime 
T 3> J a similar interplay of the Gaussian peak in the 
spectral function with the minimum of the prefactor was 
discussed in Ref. [26l . 
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APPENDIX A: ANTICOMMUTATION OF 
OPERATORS dHJ 

Since electrons are fermions, in one dimension their 
field operators are expected to satisfy the following anti- 
commutation relations: 

tpj{x)il'Y{y) + i'Y{y)ipj(x) = 0, (Ala) 
i^l{x)^l,{y) + ijl,{y)i^l{x) = 0, (Alb) 
V'7(a;)'0y(2/) + V'y(2/)'07(a;) = S.y^'S{x - y). (Ale) 

Here we check that our form of electron operators ([8]) is 
consistent with the relations (|A1|) . 

We start by discussing the commutation relations of 
operators Zi^.y and zj ^. By definition these operators 
act on the spin chain ^ and change the number of sites 
as follows. Operator Zi^j removes site I from the spin 
chain if that site has spin 7 and gives zero otherwise. 
Conversely, the operator Zj ^ adds a new site / with spin 
7 to the spin chain by inserting it between the sites / — 1 
and 

Let us consider the effect of operator Zi^jZii^yi with 
I < I'. This operator first removes spin 7' from site I' 
and then spin 7 from site I. Alternatively, one can first 
remove site / and notice that the numbering of all sites 
after I has shifted by 1. Thus to achieve the same result, 
at the second step one needs to remove site I' ~ 1. We 
therefore conclude 

Zi^jZii^ryi = Zii^i^j/Zi^j, I < I'. (A2a) 

Repeating these arguments with operators we find 

zIzIy = 4+i^yZIv (A2b) 
Zi^^zl^, = zl_^ .^,Zi^^, I < I'. (A2c) 

Let us check Eq. (|Alc|) at x < y. We start by writing 
the first term as 

M^)^l'iy) - ^i^)Zii.)nZly),y^Hy) (A3) 

and noticing that according to the definition Q the holon 
created at point y is counted in l{y), but not in l[x). Thus 
l{x) < l{y) and one can use Eq. (|A2cp . 

Using the definition ((T]) again, we find that at a; < 
y the operators ^/(a;),7 and ^'^(y) commute, while 
*(x)Zj^)_^ = Zl^^y^ ,^,■^{x). Thus we conclude 

i^i{x)i^\'{y) = 4^),^,*(x)vE't(y)Z,(,),^ 

= -zly^n''^\y)^{x)Zi^.),, 

in agreement with Eq. (jAlcp . One can easily perform a 
similar check of Eqs. (jAlap . (jAlbp . and the case of a; > y. 
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At X = y the relations (jAlaP and (jAlbp for the op- 
erators ([5]) are trivially satisfied, because '^{x)'^{x) — 
'i\x)^\x) = 0. On the other hand, the relation (jAlcP 
is less straightforward. From Eq. (|A3[) at x = y we get 



(A4) 



il}^{x)il}i,{x) = ^{x)¥{x)5^^,, 



because the operator Zi^^ZI^, first creates a site with 
spin 7' and then removes the same site with spin 7. The 
second term in the left-hand side of Eq. (jAlcP becomes 



^l,{x)^^{x) = ^\x)^{x)ZJ^^^^^,Z,^^^^^. 



(A5) 



Contrary to the expectation based upon Eq. (jAlcp . 

The reason for this apparent discrepancy is that our 
operators ^ act in a restricted Hilbert space, where two 
electrons cannot occupy the same point a;, even if their 
spins are opposite. This is a fundamental feature of our 
theory, which reflects the fact that electrons repel each 
other very strongly. In this restricted space, the operator 



'4^7{x)ipl'{x) +iIjI,{x)iP^{x) 



(A6) 



is not equivalent to 5^^/. [Here it is convenient to view 
X as Qb discrete coordinate and replace S{x — x) — > 1 in 
Eq. (|XTc)) ]. Indeed, in our Hilbert space the state at point 
X can be either empty, or occupied with a single electron, 
with possible spins t or |. When acting on these states, 
the operator (jA6[) has the following effect: 



|0) 

IT) 
li) 



"77 ' 



It is easy to check that the sum of operators (|A4[) and 
()A5|) has exactly the same effect on these three states. 
Thus our expressions ([5]) for the electron operators have 
correct anticommutation relations. 



APPENDIX B: BOSONIZATION OF THE 
OPERATORS Zi,^ IN THE CASE OF 
NON- VANISHING MAGNETIZATION 



In Sec. nil Bl we bosonized the operators Z; in the 
SU(2)-symmetric case, when no magnetic field is applied 
to the system. In the presence of the field, the ground 
state of the system has unequal densities of electrons 
with spins t and |, giving rise to a finite magnetization 
m = 2(5f). As a result the Fermi sea of the Jordan- 
Wigner fermions (Fig. ^ expands to accommodate their 
increased density (ajoj) = {1 + m)/2, see Eq. The 
Fermi points corresponding to this density are 



-(3 + m). 



(Bl) 



Our bosonization procedure in Sec. IIII Bl was performed 
for m = 0, but the derivation can be easily generalized 



to the case of m > by using the proper values (|B1|) 
of qL,B. instead of Eq. (P^ . The resulting bosonization 
expressions for the operators Z;.^ are given by 



e»f (lT"i)igTi[¥'(0 + 5(l=F"i)'?(0] 



with the upper and lower signs corresponding to 7 = "f 
and J,, respectively. 

It is instructive to substitute Eq. (jB2p into our expres- 
sion (|14p for the electron annihilation operators and com- 
pare the resulting bosonization formulas with the stan- 
dard expressions As we discussed in Sec. IIII C[ 
the right-moving electron is constructed out of a right- 
moving holon and a left-moving spinon. We therefore 
combine the second term in Eq. (|B2|) with the first term 
in the parentheses in Eq. (fT4|) and obtain 



ie(x) 



-[kpx+(t,(x)] ±i[ip(nx)-i^i)(nx)] 



v27ra 



(B3) 

This result should be compared with the standard 
bosonization expression (|15ap . 

As expected, instead of kp — kp/2 we find that the 
Fermi momentum is now a function of the magnetization: 



(B4) 



In addition, by comparing Eq. (jB3p and the analogous 
expression for the left-moving electron 



-ie{x) 



il^[k'^x+ct>(x)]^^i[ip{nx) + i^-d(nx)] 



(B5) 

with Eq. (fT5|) one finds the following relations between 
the bosonic fields: 



cbix) 

e{x) 

m 



\/2 0p(x), 
9p[x) + ni9c[x) 

4>a{l/n) ~ m4)p{l/n) 
ail In). 



(B6) 
(B7) 

(B8) 
(B9) 



These relations generalize our earlier expressions (|18b ..b) 
and (|29h .b) to the case of non-vanishing magnetization. 
It is worth noting that at m > the original charge and 
spin boson modes are mixed. This mixing was discussed 
in Ref . i51i . 



APPENDIX C: SOLUTION OF EQUATION ([Ff]) 
WITH CONDITIONS (|BT|) 

In this Appendix we show that the solution of Eq. ([ST]) 
with conditions ([CTjl has the form of Eqs. (|62ap and (|62bp . 
We first rewrite the integral equation (|57p in terms of 
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complex variables w = e*"^ and z = e^'^ . Using the fact 
that cj^ is even function of q, we find 



c (w) 



+ 1 
2w 



I , , w f dz — 1 I 
2m z — w z^ 



(CI) 

Here we use the notations c^iji) — c^(vS)\ the integral 
is taken over the unit circle \z\ = 1 in counterclockwise 
direction. 

According to the condition (j61bp function c~ (w) van- 
ishes when w is on the right semicircle R (defined as 
w = e*"? with q between —tt/2 and tt/2). In additions, 
c+ (z) vanishes on the left semicircle, see Eq. (I61ap . Thus 
for w R we have 



+ 1 



c+(w) = — [ 

Jr 



dz — 1 



z — w z^ 



c+(z), w ^ R. 



Here we assume that the contour R is traversed in the 
counterclockwise direction, from z — ~i to z — i. 

This equation can be solved using the theory of singu- 
lar integral equations42. It will be convenient to introduce 
a new unknown function 



z^ - 1 
(z)^^c+(z). 



(C3) 



and does not diverge faster than l/{w zL i) at the ends 
of the contour, is unique up to an arbitrary numerical 
coefficient. It can be found using the techniques discussed 
in S79 of Ref. S. The solution is 



1 / 1 /■ Inz , 
■ exp — / dz 

\Tn JjiZ~w 



1 



(CIO) 



where the logarithm is defined with the branch cut along 
the negative real axis. It is easy to check directly that 
the function (jClOp does satisfy the above conditions. 
For w e R, by combining Eqs. JCH, (EZl), and (IC9l) 



one finds c^{w) = $+(w). Then, substituting w = e*^ 
one finds the result (|62ap . 

Our next goal is to find c~(w) for w on the left unit 
semicircle L, defined by w = e** with q between 7r/2 and 
3tt/2. Using Eq. (|C1[) and noticing that c+(w) = on L, 
we find 



27ri ./fl z — w z^ 



■ c^(z) = —w^{w), 



(Cll) 

see Eqs. ((C3t and (IC5t. Then, substituting w = e"?, 
we obtain the result (|62b|) with the same normalization 
constant as in Eq. (|62ap . 



Then the integral equation (jC2[) takes the form 



+ 1 
w"^ — 1 



1 



(w) — — : dz 



(^) 



z — w 



weR. (C4) 



In solving this equation we will assume that the unknown 
function (t>{w) is analytic with possible exception of an 
integrable singularities at the ends z = ±i of the contour 
R. Let us introduce a new function 



^(w) = — / dz 



z — w 



(C5) 



Obviously $(«;) is analytic everywhere except the con- 
tour R and approaches zero at infinity. At the contour 
$(w) has a branch cut. When w approaches R from the 
left or right, the function $(w) takes the values 



$^(w) = — / dz ■ 
2Tn 



6 +0. (C6) 



ifi z — w ± (5 ' 
One can easily see that for w e i? we have 

$+(u)) - $"(u)) = (t){w), (C7) 
^+{w)+^-{w) - — / dz^^. (C8) 

TTlJji Z-W 



Substituting these relations in Eq. (jC4|) . we find that the 
values of (j){w) on the two sides of the branch cut satisfy 
the following linear relation: 



<I>+(w)=w^<i> (w), w e R. 



(C9) 



An analytic function that satisfies these conditions on the 
two sides of the contour R, falls off to zero at infinity. 



APPENDIX D: MAPPING OF THE BOSONIZED 
HEISENBERG SPIN CHAIN TO THE SPIN 
SECTOR OF WEAKLY INTERACTING 
ELECTRON SYSTEM 

1. Consequences of the spin-rotation symmetry 



In Sec. IIV C 31 we have evaluated the first correc- 
tion ([75)1 to the correlator c}j{q) above the Fermi point 
q — 7r/2, where the simple bosonization result (|67ap van- 
ishes. The correction originated from two sources. First, 
we found the contribution due to the deviation of the 
quadratic part of the Hamiltonian from the fixed point, 
JC 7^ 1/2, and then we included the sine-Gordon term 
(I28bp . The latter correction turned out to be larger than 
the former one by a factor of two. Here we show that this 
is a result of the spin-rotation symmetry of the problem. 

To this end we utilize the equivalence of the bosonized 
Hamiltonians of the Heisenberg spin chain ((28)) and the 
spin part (|16bp of the Hamiltonian of weakly interacting 
electrons. The exact form of the electron-electron inter- 
actions does not affect the general form of the Hamilto- 
nian (|16p . In the simplest case, one can consider only 
backscattering of electrons by each other. For the elec- 
trons in the vicinity of the Fermi level the most general 
form of backscattering Hamiltonian is 



with the coupling constant 

g^l'SS' — 7:5l(^77"5<5<5' + ""77' ■ 



(Dl) 



(D2) 
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where ct^^i is the vector of standard PauU matrices and 
gi is the 2A;p- Fourier component of the interaction po- 
tential. In Eq. (jDl[l the summation over repeating spin 
indices is imphcd. 

The exact form of the coupling of electron spins in the 
Hamiltonian ()Dip . (ID2|) is dictated by the SU(2) symme- 
try of the problem with respect to the rotation of electron 
spins. Ignoring this symmetry for the moment, we will 
view Eq. (|D2p as a special case of tensor 



'SS' 



'SS' 



(D3) 



Unlike Eq. (|D2[) . this form of coupling violates the spin- 
rotation symmetry, unless gix = giy — giz- The standard 
treatments of weakly-interacting electrons systems, in- 
cluding the derivation of the bosonizcd Hamiltonian (|16p . 
start with two constants, gi\\ = giz and gi± — gix — giy, 
and eventually equate g^ and gi±. 

Interaction constants gip, gi\\ and gi± affect different 
terms of the bosonized Hamiltonian (jl6p . Parameter gip 
corresponds to density-density coupling and affects the 
Hamiltonian of the charge degrees of freedom via renor- 
malization of Vp and Kp. It does not affect and for 
our purposes can be ignored. The coupling constant gi\\ 
enters via Kcr = 1 -f .gi||/2, whereas the spin flip scatter- 
ing accounted for by gi± transforms to the sine-Gordon 
term in Eq. (jl6bp . 

Our calculation of the correction ((73)) neglected the 
sine-Gordon term. Thus g^ in Eq. (|73p is in fact gf^. 
Since c^{q) is invariant with respect to spin rotations, 
identical contributions should be expected from coupling 
constants gix and giy. While evaluating the correction to 
c+(g) due to the sine-Gordon term (|28bp . we accounted 
for both gix and giy and, as expected, obtained twice the 
resuh (1731). 



2. Alternative evaluation of the correlators c^{q) 
near the Fermi points 



The mapping of the bosonized Hamiltonian ((28|) of 
the Heisenberg spin chain and the Hamiltonian Ha^, 
Eq. (|16bp . describing the dynamics of the spin sector of 
weakly interacting electron gas, enables one to obtain 
an alternative expression for the correlators c^{q) near 
the Fermi point. Let us consider the Green's function 
of right-moving electrons Gii{x,t) traced over the spin 
indices. It is well knownpi^ that asymptotically at large 
X and t it separates into a product of charge and spin 
factors, 



GR{x,t) = -gp{x,t)gaix,t). 

TT 



(D4) 



The two factors are most easily computed using the 
bosonization transformation (jlSp . in which case the 
charge and spin factors are obtained by averaging the ex- 
ponentials of bosonic fields 0p, 9p and ^o-, 0^ in Eq. (|15ap . 




FIG. 4: The second-order diagram for the perturbative cal- 
culation of the electronic Green's function GR{x,t). 



respectively. In the absence of electron-electron interac- 
tions the parameters of the Hamiltonian take unper- 
turbed values Vp = Va- = vp, Kp = K^r = 1, and gi±^ — 0. 
Then one finds 



gf\x,t) ^ g^^\x,t) = 



1 



{x — vpt + iSsgnty^^ ' 



(D5) 



and Eq. (jD4p reproduces the standard expression for the 
Green's function of non-interacting electrons. 

In the presence of interactions the parameters of the 
Hamiltonian (jl6p renormalize, and the charge and spin 
components of the Green's function show non-trivial be- 
havior. We showed in Sec. IIII Cl that the exponentials of 
the bosonic fields (jj^ and 9a- in Eq. (jl5ap are equivalent 
to the bosonized expression (|4T|l for the operators Zi^j. 
Thus the correlator {Zi^^Z\^ in the definition (|51ap of 
c^isi) can be found from ga{x,t) at x = l/n and t — 0. 

We will find g„{x,t) by calculating electron Green's 
function (|D4p and identifying its spin component. Since 
the Hamiltonian (jl6bp of the spin sector is universal at 
low energies, the specific form of the electron-electron 
interactions is not important. It is most convenient to 
choose the form (|D1I) with the coupling constant (|D3p 
chosen so that gip = and gix = giy = giz = gi- In this 
case the interactions do not affect the charge sector, and 
gp{x,t) retains its unperturbed value (|D5p . and 



ga{x, t) — 'k{x — Vpt -I- i5 sgnt)^ ^"^ G pi{x , t). 



(D6) 



The electronic Green's function can be studied using 
straightforward perturbation theory in the coupling con- 
stant gi. 

The first non-vanishing contribution to Gpi{x^t) ap- 
pears in the second order and is given by the diagram 
Fig. [4] The respective correction to the Green's function 
is 

5GR{k,uj) = 2Gf{k,uj)J:{k,uj)Gf(k,uj), (D7) 

where the factor of 2 accounts for the trace over spin 
variables and the self-energy 

v^i, ^ ^dl f dkidk2duJiduj2 {0),, ^r^(o)/, n 
^[k,uj) ^ — J (2^;;^4 (fc2,t^2) 

xGf{k + ki~k2,UJ + UJl-UJ2). (D8) 
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Here the factor Sgf appears as the sum gf^ + g\y + g\^ 
Substituting the unperturbed Green's functions 



1 



(D9a) 



uo + vpk — iS sgnfc 
Gf{k,uj) = -, (D9b) 

and performing integration with respect to wi, [^2, and 
one of the momenta, we find 



I](fc,u;) = S„(fc) + S(fc,u;), 



where 



(DIG) 



(DU) 



B{k^q) 



uj — vpk — 2vpq + iS 
°"^dq. (D12) 



iglito — Vpk) 
IQn'^vp jQ 
ej-k + q) 
LU — Vpk + 2vpq — iS 



Here we have introduced the short-distance cutoff a for 
the electron-electron interactions. In principle, this cut- 
off may not coincide with the bandwidth cutoff a used 
in the bosonization procedure. This distinction is not 
important for the present discussion. 

The two second-order contributions to the electron self- 
energy have very different meanings. The term Sy ac- 
counts for a small correction to the velocity of spin exci- 
tations, which for our purposes can be ignored. On the 
other hand, S leads to the logarithmic renormalization of 
the electron Green's function, which affects the singular 
behavior of the correlators c^{q) near the Fermi points. 
We therefore explore this correction in more detail. 



a. Logarithmic correction to the Green's function 



Let us now substitute the expression (jD12|l for the self- 
energy I](fc, uj) in Eq. (|D7[) and perform the Fourier trans- 
formation to X and t variables. The resulting correction 
to the Green's function has the form 



SGR{x,t) = 



35 



1 



32n^Vp X — vpt + iS sgnt 



X In ■ 



{a + ivp\t\)'^ + x'^ 



(D13) 



This expression is consistent with the logarithmic renor- 
malization of the electron Green's function studied ear- 
lier in the g-ology theory, cf. Eq. (4.24) of Ref. 0. (Our 
Eq. (|D13[) is obtained by neglecting the constant g2 and 
replacing 4gf with Sgf to account for the fact that we 
assume gip — 0.) 

We can now separate the spin component g„{x,t) of 
the electron Green's function Gu{x, t) following the pre- 



scription (|D6p and obtain 



{x — Vpt + i5 sgni)i/2 

Zyl {a + ivp\t\f + x"^ 



X 1 



32 



In- 



(D14) 



where yi — gi/irvp. As expected, gs{x,+0) reproduces 
the logarithmic correction to the correlator {ZijZq ^) 
used in Eq. ([72]) . with the additional factor of 3 correctly 
included. 



6. Renormalization of the Green's function 



The logarithmic correction (jD14p to the Green's func- 
tion grows at long distances, and can in principle be- 
come large despite the smallness of the prefactor yf. To 
find out whether this is the case, one can compute the 
Green's function in the leading logarithm approximation. 
We accomplish this by adopting the multiplicative renor- 
malization procedure of Ref. I44j. We present the spin 
component of the Green's function as 

g„{x,t) = d{Og^°Hx,t). (D15) 

Here the new function c?(^) is expected to depend on x 
and t very slowly, via their logarithm. 



2 



ivp\t\y 



(DI6) 



In particular, our result (ID14P has the form (jD15|) with 



1 



16 



(D17) 



As we increase ^, the correction (ID17P grows and may no 
longer remain small. In addition, the coupling constant 
yi itself depends on ^ as 

yi 



(DI8) 



see Ref. Hi 

Following the general prescriptionii to account for the 
multiplicative corrections to the Green's function, we use 
Eq. ()D17|) to write the renormalization group equation 
upon Ind, 



lnd(0 
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(D19) 



Solution of this equation with yi(^) given by Eq. (|D18p 
has the form 



d{0 = exp - — 



yU 

16 1 -I- yiC 



(D20) 



At small yiS, it reproduces the perturbative expansion 
(jD17p , while in the limit ^ — > cx) we obtain a finite renor- 
malization of the Green's function. 



ga{x,t) = exp — — 



16 



1 



{x — Vpt + «(5sgnt)i/2 



(D21) 
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In a weakly interacting electron gas the coupling con- 
stant 2/1 <C 1, and the renormalization (|D21|) can be ig- 
nored. As interactions become stronger, yi increases and 
reaches values of order unity. Thus one can expect that 



the coefficient x in the asymptotes (|75p of the correla- 
tors (q) near the Fermi point will slowly decrease from 
its numerically obtained value x ^ 0.8 to a significantly 
lower number at q — 7r/2. 
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